Antimatter production in central Au+Au collisions at v /snn=200 GeV 
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We have used the dynamically constrained phase space coalescence model to investigate the pro- 
duction of light nuclei (anti-nuclei) based on the 1.134xl0 7 hadronic final states generated by the 
PACIAE model for the 0-5% most central Au+Au collisions at -^/snn=200 GeV with \y\ <1 and 
Pt <5 acceptances. The STAR data of !Lff, aH, 3 He, and 3 He yields and ratios are well reproduced 
by the corresponding PACIAE results. The transverse momentum distribution of z He ( 3 He) and 
j-H (\H) is also given. It turned out that the transverse momentum distribution of light nuclei is 
close to that of the corresponding anti-nuclei. 



PACS numbers: 25.75.-q, 24.85. +p, 24.10.Lx 

I. INTRODUCTION 

The nucleus- nucleus collisions at top RHIC energy 
produce an initial hot and dense matter (quark-gluon 
mater, QGM). It has been interpreted as a strongly cou- 
pled quark-gluon plasma (sQGP) [H-Q. This is nearly a 
perfect liquid composed of quarks and gluons but is not 
a free gas-like quark-gluon plasma (fQGP) expected by 
theorists and experimentalists long time ago. 

The anti-nuclei production is great importance in the 
nuclear and particle physics, the astrophysics, and the 
cosmology. One believes that the matter and antimatter 
exist in equal abundance during the initial stage of the 
universe. However, a mystery exists: how this symmetry 
got lost in the evolution of the universe with no signif- 
icant amount of antimatter being present. Because the 
initial fireball created in ultra-relativistic heavy ion col- 
lisions is similar to the initial stage of the universe, the 
study of anti-nuclei production in ultra-relativistic heavy 
ion collisions may light this issue. However, the study of 
light nuclei (anti-nuclei) production is quite hard both 
experimentally and theoretically because of the low pro- 
duction multiplicity. 

The STAR collaboration has reported their measure- 
ments of \H and j-H in Au+Au collisions at the top 
RHIC energy [5J. They have measured "70+17 antihy- 
pertritons (^H) and 157+30 hypertritons (\H)" in the 
89 million minimum-bias and 22 million central ( "head- 
on") Au+Au collision events at y's NN =200 GeV. Thus 
the corresponding yields are estimated to be 6.31 xlO -7 
and 1.41xl0~ 6 , respectively. The ALICE collaboration 
has also published their preliminary d yield of~6xl0~ 5 
measured in the pp collisions at TeV 0,0|. 

On the other hand, the theoretical study of light nu- 
clei (anti-nuclei) is usually separated into two steps. The 
nucleons and hyperons are first calculated with some 
selected models, such as the transport models. Then 
the light nuclei (anti-nuclei) are calculated by the phase 
space coalescence model |8l-tlfJj and/or the statistical 



model [U EH etc. Recently, production of light nuclei 
(hypernuclei) in Au+Au/Pb+Pb collisions at relativis- 
tic energies have been investigated theoretically by the 
coalescence+blast-wave method [l3[ and the UrQMD- 
hydro hybrid model+ thermal model [l4j, respectively. 

We have proposed an approach studying the light nu- 
clei (anti-nuclei) production in ultra-relativistic pp col- 
lisions by dynamically constrained phase-space coales- 
cence model [15J. This approach is based on the final 
hadronic state generated by a parton and hadron cascade 
model PACIAE The calculated light nuclei (anti- 

nuclei) yield in non-single diffractive (NSD) pp collisions 
at y/s=7 TeV is well comparing with ALICE data [1,[3|- 
The transverse momentum distribution and rapidity dis- 
tribution are also predicted for 3 He ( 3 He) and^-if 

in NSD pp collisions at y / s=7 and 14 TeV. In this paper, 
we use this method to investigate the light nuclei (anti- 
nuclei) and hypernuclei (anti-hypernuclei) productions 
in the Au+Au collisions at ^/snn =200 GeV. 

The paper is organized as follows: In the Sec. II, we 
briefly introduce the PACIAE model and the dynami- 
cally constrained coalescence model. In Sec. Ill, the cal- 
culated light nuclei (anti-nuclei) and hypernuclei (anti- 
hypernuclei) yields, ratios, as well as the transverse mo- 
mentum distributions are given and compared with the 
STAR data. A short summary is the content of Sec. IV. 



II. MODELS 

The PYTHIA model [l7| is devised for the high en- 
ergy hadron-hadron (hh) collisions. In this model, a 
hh collision is decomposed into the parton-parton colli- 
sions. The hard parton-parton scattering is described by 
the leading order perturbative QCD (LO-pQCD) parton- 
parton interactions with the modification of parton dis- 
tribution function in a hadron. The soft parton-parton 
collision, a non-perturbative phenomenon, is considered 
empirically. The initial- and final-state QCD radiations 
and the multiparton interactions are also taken into ac- 
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count. Therefore, the consequence of a hh collision is 
a partonic multijet state composed of di-quarks (anti- 
diquarks), quarks (antiquarks) and gluons, as well as a 
few hadronic remnants. This is then followed by the 
string construction and fragmentation. A hadronic final 
state is obtained for a hh collision eventually. 

The parton and hadron cascade model PACIAE [l|| is 
based on PYTHIA and is devised for the nucleus-nucleus 
collisions mainly. In the PACIAE model, first of all the 
nucleus-nucleus collision is decomposed into the nucleon- 
nucleon (NN) collisions according to the collision geom- 
etry and NN total cross section. Each NN collision is 
described by the PYTHIA model with the string frag- 
mentation switches-off and the di-quarks (anti-diquarks) 
randomly breaks into quarks (anti-quarks). So the con- 
sequence of a NN collision is now a partonic initial state 
composed of quarks, anti-quarks, and gluons. Provided 
all NN collisions are exhausted, one obtains a partonic 
initial state for a nucleus-nucleus collision. This par- 
tonic initial state is regarded as the quark-gluon matter 
(QGM) formed in the relativistic nucleus-nucleus col- 
lisions. Secondary, the parton rescattering proceeds. 
The rescattering among partons in QGM is randomly 
considered by the 2—^2 LO-pQCD parton-parton in- 
teraction cross sections [lq . In addition, a K factor 
is introduced here to include the higher order and the 
non-perturbative corrections. Thirdly, the hadroniza- 
tion follows after the parton rescattering. The partonic 
matter can be hadronized by the Lund string fragmen- 
tation regime |l7ll and/or the phenomenological coales- 
cence model [la]. Finally, the hadronic matter proceeds 
rescattering until the hadronic freeze-out (the exhaus- 
tion of the hadron-hadron collision pairs). We refer to 
[l6j for the details. 

In quantum statistical mechanics [l9| one can not pre- 
cisely define both position q = (x,y,z) and momentum 
p = ipxiPyiVz) of a particle in the six dimension phase 
space, because of the uncertainty principle 

AqAp~ h 3 . 

We can only say this particle lies somewhere within a 
six dimension quantum "box" or "state" with volume 
of AqAp. A particle state occupies a volume of h 3 in 
the six dimension phase space [19(. Therefore one can 
estimate the yield of a single particle by 
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1 if 1 = p, 2 = fi, 3 = A; 
j m ^ m inv ^ to + Am; 

~ s |gia|<A>, l&sKA), |&»|<A>; 

otherwise, 

(4) 

where 

m inv = [(Eh. +E 2 + E 3 f - (p a + p 2 + p 3 ) 2 } 1/2 , (5) 

here (Ei, E2,E 3 ) and (pi,p2,ps) are the energy and mo- 
mentum of particles p, n, A, respectively. In Eq. (j4]), too 
and Do stand for, respectively, the rest mass and diam- 
eter of j-H, Am refers to the allowed mass uncertainty, 
and \qij\ = \<jfc — <jj\ is the vector distance between parti- 
cles i and j. 

III. CALCULATIONS AND RESULTS 

As the hadron position and momentum distributions 
from transport model simulation are discrete, the inte- 
gral over continuous distributions in Eq. (j4|) should be 
replaced by the sum over discrete distributions. In a 
single event of the final hadronic state obtained from 
transport model simulation, the configuration of 
(p+n+A) system can be expressed as 

Cp n A(Qi,Q2,q3;Pi,P2,P3), 



(6) 



where the subscripts 1 = p, 2 = n, 3 = A, and q\ refers 
to the distance between p and the center-of-mass of p, n, 
and A for instance. Then the third constraint (diameter 
constraint) in Eq. Q is correspondingly replaced by 

qi < Ro, q 2 < Ro, <73 < Ro, (7) 

where Rn refers to the radius of ^-H. 

_ A 

Each of the above configuration contributes a partial 
yield of 



(1) 



1 if too ^ irii nv ^ too + Am, 
'in-, = \ qi < Ro, q_2 < Ro, 93 < -Ro; 

otherwise; 



(8) 



where H and E are the Hamiltonian and energy of the 
particle, respectively. Similarly, the yield of N particle 
cluster can be estimated by 



Y N - 



dq\dpi...dq^fdp 



N 



H<E 



h 3N 



(2) 



Therefore the j-H yield in our dynamically constrained 
phase space coalescence model, for instance, is assumed 



to the j^H. So the total yield of ^-H in a single event 
is the sum of the above partial yield over the configura- 
tions of Eq. ([6]) and their combinations. An average over 
events is required at the end. 

In the PACIAE simulations we assume that the hy- 
perons heavier than A already decay. The model param- 
eters are fixed on the default values given in PYTHIA. 
However, the K factor as well as the parameters parj(l), 
parj(2), and parj(3), relevant to the strange production 
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TABLE I: Strange particle rapidity density dN/dy at the 
mid-rapidity (\y\ <1 for A and A, \y\ <0.75 for H~ and H~) in 
the 0-5% most central Au+ Au central collisions at % /snn=200 
GeV. 



Particle type 


STAR" 


PACIAE 


A 


16.7 ± 0.2 ± 1.1 


16.8 


A 


12.7±0.2±0.9 


13.1 




2.17±0.06±0.19 


2.08 




1.83±0.05 ± 0.20 


1.68 



a The STAR data were taken from [20] 



in PYTHIA [T^, are roughly fitted to the STAR data of 

A, A, S - , and 5 in the 0-5% most central Au+Au col- 
lisions at v / sii7=200 GeV as shown in Tab. U The 
fitted parameters of K=3 (default value is 1 or 1.5 fjl\). 
parj(l)=0.12 (0.1), parj(2)=0.55 (0.3), and parj(3)=0.65 
(0.4) are used to generate 1.134x 10 7 final hadronic states 
by the PACIAE model for the 0-5% most central Au+Au 
collisions at ^/snn =200 with \y\ <1 and p t <5 accep- 
tances. Then d (d), 3 He ( 3 Se), as well as (pf) 
yields and ratios are calculated by the dynamically con- 
strained phase-space coalescence model. The calculated 
results are denoted as "PACIAE" later. 

Figure Q] shows the calculated hadron yields (solid tri- 
angles) in the Au+Au collisions at */snn=200 GeV. The 
open circles in this figure are the experimental data taken 
from [e|. One sees in this figure that the PACIAE re- 
sults well agree with the experimental data. 
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FIG. 1: Hadron yields (including d,d) in the Au+Au collisions 
at ^/snn=200 GeV. The open circles are the experimental 
data taken from pj| and the solid triangles are the PACIAE 
results. 



In the Tab. HQ d, d, 3 He, 3 He, 3 A H, and \H yields 
are given. The STAR yield of 3 He ( 3 He) is estimated 
by the \H (pf) yield of 1.41E-06 (6.31E-07) divides 
by the ratio of \H to 3 He (pf to \H) of 0.82 (0.89) 



TABLE II: Light nuclei (anti-nuclei) rapidity density dN/dy 
at the midrapidity (\y\ <1) in the Au+Au collisions at 
^5^=200 GeV. 

Nucleus Exp, data PACIAE Ref. [13]° Rcf. [14] b 

d 0.098 c 0.085 d 

d 0.045 c 0.041 d 

\H 1.41E-06 e 1.15E-06 / >1.05E-04 4E-05 

OF 6.31E-07 e 5.29E-07 / >4.90E-05 

A 

3 He 1.72E-06 9 l.SOE-OG-f >1.65E-04 

3 He 7.09E-07' 1 6.17E-07 / >7.30E-05 
a taken from Tab. I in [13]. 
6 taken from Fig. 6 in [14]. 
c taken from [121 . 

d calculated with Am=0.0003 GeV. 

e estimated from [5J. 

f calculated with Am=0.00015 GeV. 

9 equal to 1.41E-06/0.82 (0.82 is taken from [5]). 

h equal to 6.31E-07/0.89 (0.89 is taken from [5]). 

0. We see in this table that the agreements between 
experimental data and PACIAE results are well. 

We give the ratios of anti-nuclei to nuclei in Tab. IIIII 
One sees again in this table that the PACIAE results are 
well comparing with the STAR data. 

TABLE III: Light anti-nuclei to nuclei ratios in the Au+Au 
collisions at v/ s nn=200 GeV. The STAR data are taken from 
1- 



Ratio 


STAR 


PACIAE 


3 ~He/ 3 He 


0.45+0.18 + 0.07 


0.41 




0.49+0.18 + 0.07 


0.46 


\H/ 3 He 


0.82+0.16 + 0.32 


0.76 


A ' 


0.89+0.28 + 0.13 


0.86 



We plot in Fig. [2] the calculated d, d, \H, LH, 

3 He, and 3 He transverse momentum distributions in 
the Au+Au collisions at ^/snn=200 GeV. Figure 0(a), 
(b), and (c) are the calculated p t ditribution for d (d), 
{j-H), and 3 He { 3 He), respectively. The pattern 

of 3 He transverse momentum distribution is consistent 
with the pt distribution of 3 He in the Pb+Pb collisions 
at y/sNN =2. 76 TeV 0- The peak of p t distribution in 
the Pb+Pb collision is located at higher p t than the one 
in the Au+Au collision because of the different reaction 
energy. The strong fluctuation showing in the panel (b) 
and (c) indicates that the 1.134xl0 7 events are still not 
enough for these pt distributions. 

The PACIAE results of light nuclei (anti-nuclei) av- 
erage transverse momentum in the Au+Au collisions 
at v /snn=200 GeV are given in Tab. [TV] Here we see 
that the average transverse momentum of light nuclei 
is nearly equal to the one of corresponding anti-nuclei. 
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FIG. 2: (Color online) PACIAE results of the transverse momentum distributions of light nuclei (anti-nuclei) in the Au+Au 
collisions at ^/snn=200 GeV. In this figure the blue dashed histograms are calculated for nuclei and red solid histograms for 

anti-nuclei. The panel (a) is calculated for d and d, (b) for ^-H and \H, and (c) for 3 He and 3 He. 



TABLE IV: PACIAE results of the light nuclei (anti-nuclei) 
average transverse momentum (p t ) in the Au+Au collisions 
at yi^=200 GeV. 



D 


D 


3 He 3 He \H 




0.92 


0.96 


1.05 1.06 1.15 


1.18 



This feature is already seen in the hadron (anti-hadron) 
production. 

IV. CONCLUSION 

In summary, we have employed the dynamically con- 
strained phase space coalescence model to investigate the 
light nuclei (anti-nuclei) production based on the final 
hadronic state generated by the PACIAE model for the 
0-5% most central Au+Au collisions at v /snn=200 GeV 
with \y\ <1 and <5 acceptances. The calculated PA- 



CIAE results of d (d) yield of 0.085 (0.041) is close to the 
experimental datum of 0.098 (0.045) 0. The PACIAE 
yields of j-H, \H, 3 He, and 3 He as well as their ratios 
are also consistent with the STAR data [5] . The PACIAE 
results of light nuclei (anti-nuclei) transverse momentum 
distributions are also given. The consistency between 
the PACIAE results and the corresponding experimen- 
tal data demonstrates that the PACIAE+dynamically 
constrained phase space coalescence method is able to 
describe the production of light nuclei (anti-nuclei) and 
hypernuclei (anti-hypernuclei) in the relativity heavy ion 
collisions. 
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